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^ ABSTRACT 

The magnetic field of the classical T Tauri star V2129 Oph can be modeled ap- 
^ proximately by superposing slightly tilted dipole and octup ole moments, wi th polar 

• magnetic field strengths of 0.35 kG and 1.2 kG respectively ( Donati et al. f2007, here- 

^ after D07). Here we construct a numerical model of V2129 Oph incorporating this 

result and simulate accretion onto the star using a three-dimensional magnetohydro- 
Q dynamic code. Simulations show that the disk is truncated by the dipole component 

and matter flows towards the star in two funnel streams. Closer to the star, the flow 
c/^ is redirected by the octupolar component, with some of the matter flowing towards 

, the high-latitude poles, and the rest into the octupolar belts. The shape and position 

of the spots differ from those in a pure dipole case, where crescent-shaped spots are 
^ observed at the intermediate latitudes. 

^ Simulations show that if the disk is truncated at the distance of r ^ 6.2R^ which 

1— I is comparable with the co-rotation radius, rcor ^ Q.SR^, then the high-latitude polar 

OO spot dominates, but the accretion rate obtained from the simulations (and from the 

accompanying theoretical calculations) is about an order of magnitude lower than 
^ the observed one. The accretion rate matches the observed one if the disk is dis- 

^ rupted much closer to the star, at 3AR^. However, in that case the octupolar belt spots 

1— I strongly dominate. In the intermediate case of r ^ 4.3i?^, the polar spots are suffi- 

ciently bright, and the accretion rate is within the error bar of the observed accretion 
O rate, and this model can explain the observations. However, an even better match has 

^ been obtained in experiments with a dipole field twice as strong compared with one 

• ^-j suggested by D07. 

/\ The torque on the star from the disk-magnetosphere interaction is small, and the 

H time-scale of spin evolution, 2 x 10^ — 10^ years is longer than the 2 x 10^ years 

age of V2129 Oph. This means that V2129 Oph probably lost most of its angular 
momentum in the early stages of its evolution, possibly, during the stage when it was 
fully convective, and had a stronger magnetic field. 

The external magnetic flux of the star is strongly influenced by the disk: the field 
lines connecting the disk and the star inflate and form magnetic towers above and 
below the disk. The potential (vacuum) approximation is still valid inside the Alfven 
(magnetospheric) surface where the magnetic stress dominates over the matter stress. 
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1 INTRODUCTION 



Classical T Tauri stars (CTTSs) are young, low-mass stars and 
often show signs of a strong magnetic field (e.g./Basri, Marcy 
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& Valenti|1992||Johns-Krull, Valenti & Kor esko 1999| [Johns' 



Krull|2007 ^. The magnetic field plays a crucial role in disc ac- 
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cretion by disrupting the inner regions of the disc and chan- 
nehng the disk matter onto the star ( [Ghosh & Lamb|[r978| 
|l979( [Konigl|[1991|. Zeeman measurements of the surface 
magnetic fields of CTTSs show that the field is strongly non- 
dipole in many stars ([Donati & Collier Cam eron||1997[ |Do-| 
Inati et aL|19"99l|Jardine et aL|2002|,20Q6;.Johns-Krull|2007P . 
The magnetic field of the CTTS V2129 Oph has been recently 
measured with the ESPaDOnS and NARVAL spectropolarime- 
ters (D07). Decomposition of the observed field into multi- 
polar components has shown that the field of V2129 is pre- 
dominantly octupolar {Boot = 1.2 kG), with a smaller dipole 
component {Bdip = 0.35 kG) and a much smaller contribu- 
tion from the higher multipole components. The large-scale 
magnetic field geometry has been calculated by extrapolating 
the surface magnetic field to larger distances by suggesting 
that the field is potential (vacuum approximation). The data 
presented in D07 correspond to the June 2005 observations. 
More recent measurements (July 2009) have shown that the 
fiel d of V2129 Oph is larger: Boct = 2.1 kG and Bdip = 0.9 
kG ( [Donati et al.|2010 , hereafter DIO). Those authors suggest 
that the field of V2129 Oph varies in time due to the dynamo 
process. In this paper we assume that the field is fixed and 
corresponds to that given by D07. However, we do perform 
test experimental runs at larger values of the dipole field. 

Analysis of the chromospheric spot distribution in V2129 
Oph have shown that the polar, high-latitude spot dominates 
(D07, see also [Jardine, Gregory & Donati] 2 008) ) . D07 devel- 
oped an approximate model where they investigated possi- 
ble paths of matter flow in the potential field along closed 
field lines (see e.g., ^Gregory et al.[[2006l ). They concluded 
that the disk should be disrupted by the dipole component 
far away from the star, at r ^ 7i?^ . In this case matter flows 
mainly along the field lines of the dipole component and hits 
the star at high latitudes which corresponds to observations. 
These estimates are useful first-order approximations to the 
processes in V2129 Oph. As a next step it is important to in- 
vestigate matter accretion onto V2129 Oph in global MHD 
simulations, where the assumption that the field is potential 
can be dropped, and the interaction of the disk with the com- 
plex field of the star can be investigated. 

In this paper, we developed a model star with parame- 
ters of V2129 Oph and with magnetic field distribution close 
to that observed in V2129 Oph (as described in D07). We ap- 
proximated the field with superposition of the main (dipole 
and octupole) moments, and performed global 3D MHD simu- 
lations of matter flow around such a star. In the past, we were 
able to model accretion onto stars with quadrupole (Long, 
Romanova & Lovelace'2007[ [2008| ), and octupole ( [Long, Ro- 
manova & Lamb 2010a) magnetic field components. In this 
paper we investigate for the first time accretion onto a star 
with realistic magnetic field and parameters corresponding to 
CTTS V2129 Oph. In the accompanying paper, we investigate 
numerically accretion onto a model star with parameters close 
to BP Tau ([ Long et al.|2010b |). 

The goals of this paper are the following: (1) to investi- 
gate magnetospheric accretion onto this model star in global 
3D MHD simulations, (2) to investigate numerically the mag- 
netic field structure around the star and to compare it with the 
initial, potential field, (3) to compare the simulated and ob- 
served accretion spots, (4) to compare the results with results 
for a purely dipolar field, (5) to obtain the accretion rate from 



simulations and compare it with that obtained from observa- 
tions, (6) to investigate the torque on the star. 

Section 2 briefly describes the numerical model we use in 
our simulations. Section 3 compares the observed and mod- 
eled magnetic fields of V2129 Oph. Section 4 discusses the 
results of the simulations. Our conclusions are summarized in 
Section 5. 



2 NUMERICAL MODEL AND REFERENCE VALUES 

We solve the 3D MHD equations with a Godunov-type code 
in a reference frame rotating with the star using the "cubed 
sphere" grid. The model has been described earlier in a series 
of previous papers (e.g., Koldoba et al. 2002{ [Romanova et al.j 



|2004a[ [Long, Romanova & Lamb|2010a| ). Hence, we describe 
it only briefly here. 

Initial conditions. A rotating magnetic star is surrounded 
by an accretion disc and a corona. The disc is cold and 
dense, while the corona is hot and rarefied, and at the ref- 
erence point (the inner edge of the disc in the disc plane), 
Tc = lOOTd, and pc = O.Olpd, where the subscripts 'd' and 'c' 
denote the disc and the corona. The disc and corona are ini- 
tially in rotational hydrodynamic equilibrium (see [Romanova| 
|et al.|[2002| for details) . An a-type viscosity is incorporated 
into the code and helps regulate the accretion rate. The vis- 
cosity is nonzero only in the disc, that is, above a certain 
threshold density (which is p = 0.2 in the dimensionless units 
discussed below). We use a = 0.04 in this work. 

Boundary conditions. At both the inner and outer bound- 
aries, most of the variables A are set to have free boundary 
conditions, dA/dr = 0. On the star (the inner boundary) the 
magnetic field is frozen into the surface of the star, that is, the 
normal component of the field, Bn, is fixed, though all three 
magnetic field components may vary. At the outer boundary, 
matter is not permitted to flow back into the region. The free 
boundary condition on the hydrodynamic variables at the stel- 
lar surface means that accreting gas can cross the surface of 
the star without creating a disturbance in the star or the flow. 
This neglects the complex physics of the interaction of the ac- 
creting gas with the star, which is expected to produce a shock 
wave in the stellar atmosphere (e.g., Koldo ba et al. 12008] ). 

A "cubed sphere" grid is used in the simulations. The grid 
consists of Nr concentric spheres, where each sphere repre- 
sents an inflated cube. Each of the six sides of the inflated 
cube has an x curvilinear Cartesian grid. The whole grid 
consists of 6 X A^r X A^^ cells. In our dipole + octupole model 
of V2129 Oph, the octupole component dominates and a high 
radial grid resolution is needed near the star. We achieve this 
by choosing the radial size of the grid cells to be 2.5 times 
smaller than the angular size in the region r < 4R^, while 
it is equal to the angular size in the outer region as in all 
our previous work. The typical grid used in simulations has 
6x NrX N'^ = 6x150 X 41^ grid cells. 

Reference values. The simulations are performed in di- 
mensionless variables A. To obtain the physical dimensional 
values A, the dimensionless values A should be multiplied by 
the corresponding reference values Ao as A = AAq. These 
reference values include: mass Mo = M^, where is the 
mass of the star; distance Ro = R^/0.3b, where R^ is the ra- 
dius of the star; velocity vo = {GM / RqY^'^ . The time scale is 
period of rotation at r = i?o : ft = 2tiRq /vq . 
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Reference values: 


M,{Mq) 


1.35 




2.4 


Ro (cm) 


4.8 X IQii 


vq (cm s~^) 


1.9 X 10'^ 


Po (days) 


1.8 



Table 1. The reference values which are fixed in all models. 




To derive the reference values for the magnetic field, we 
take into account that the tilts of the dipole and octupole mo- 
ments relative to the rotational axis are small, and hence we 
use the formulae for the axisymm etric case: /ii — O.bBi^Rl 
and fi3 = 0.25B3i.Rl (see, e.g.. Long, Romanova & Lamb 
2010ap . We suggest that the polar field: Bi^ = Bdip, = 
Boot' The reference magnetic moments for the dipole and oc- 
tupole components are /ii,o = BqRI and fi3,o = BqRq re- 
spectively, where Bo is the reference magnetic field. Hence, 
the dimensionless magnetic moments are: Jli = /ii//ii,o, 
Jl^ = /is/Ms, 0. We take one of the above relationships; for 
example, the one for the dipole; to obtain 



Bn 



Mi,o 



0.5Bi^ fR 



Ml 



Hence, the reference magnetic field depends on the dimen- 
sionless parameter Jn and the dipole magnetic field, Bi^. The 
reference density po = Bl /vl and the reference mass accre- 
tion rate Mo = pqvqRq also depend on these parameters: 



po = Bl/vl = 1.5 X 10" 



2 / \ 2 



Bi. \ ^ 
350G ) ) cm- 



Mo — poVqRq 



'Mo 
yr ■ 



(2) 



(3) 



We find the other reference values in a similar way: the torque 
No = poVqRo; energy flux, Eo = po^o^o; energy flux per 
unit area, Fq = poVq; temperature To = IZpo/po, where K 
is the gas constant; and the effective blackbody temperature 
Teff,o = {povo/ay^^^, where a is the Stefan-Boltzmann con- 
stant. Tab.[T]shows the reference values which are fixed in all 
models. Tab.[2]shows the reference values for our main mod- 
els which depend on the parameters Jli and Bi^ (see §4). In 
the paper, we show our results in dimensional units. Formulae 
|2][3]and Tab.|2]help scale the results to other parameters. 

In most of the models we change the dimensionless oc- 
tupolar moment /Is in same proportion as Jli so as to keep 
the ratio ps/pi fixed (models 1-3). To find the ratio be- 
tween the dimensionless moments, we use the approximate 
formulae for aUgned moments: and obtain the ratio ps/Jli — 
B3^Rl/2Bi^ ^ 0.22 (where R^ = 0.35). In a few test cases, 
we increased the relative strength of the dipole field (models 
4-5). 



3 RECONSTRUCTED MAGNETIC FIELD OF V2129 OPH 
AND THE FIELD USED IN 3D SIMULATIONS 

The magnetic field of V2129 Oph has been reconstructed 
from observations by D07. The authors decomposed the mea- 
sured surface magnetic field into its spherical harmonic com- 



Figure 1. Polar projections of the magnetic field components at the 
stellar surface in the dipole -I- octupole model of V2129 Oph in spher- 
ical coordinates: radial magnetic field, Br', azimuthal magnetic field, 
B^ ; and meridional magnetic field, Bq . The outer boundary, the bold 
circle and the two inner dashed circles represent the latitude of —30°, 
the equator, and the latitudes of 30° and 60°, respectively. The red 
and blue regions represent positive and negative polarities of the 
magnetic field. 



ponents. They then approximated the field as a superposition 
of low-order multipoles, adjusting the magnitudes of the mul- 
tipole moments and their orientations to give the best match 
of the modeled field to the measured field. 

In one of the best matches the surface magnetic field of 
V2129 Oph is approximated by a 0.35 kG dipole and a 1.2 kG 
octupole field, with corresponding magnetic moments tilted 
by Gi = 30° and 63 = 20° (towards phases 0.7 and 0.9) with 
respect to the rotational axis (D07).We choose this field dis- 
tribution as the one for our dipole + octupole model. Higher- 
order multipoles and azimuthal component are also present, 
but they carry much less magnetic energy than the two main 
components. 

Note that method of field reconstruction from observa- 
tions has some uncertainties, and the error bars are different 
in different stars and for different components. The error bar 
is typically ^25% for the strongest field component (the oc- 
tupole in the case of V2129 Oph) but larger for the weaker 
field components. In V2129 Oph, the dipole field is 3-4 times 
weaker than the octupole, and the error bar can be 3-4 times 
larger, that is < 100%. The error bars on the phases and tilts 
of the magnetic axes (about the rotational axis) also depend 
on the actual values: the phases are less accurate when the 
tilts are low. In the case of V2129 Oph, the error bar in the 
phase of the dipole axis is likely to be 0.1-0.2 rotation cycle, 
or larger, while the tilt of the dipole axis is accurate to no 
better than 10° - 20°. 

Fig. [1] shows the polar projection of the modeled (dipole 
plus octupole) field on the star's surface. The style of the fig- 
ures and the components of the field Br, B(p and Be are simi- 
lar to those reconstructed from observations by D07 (see their 
Figs. 12 and 13). Comparisons of our plots with the D07 plots 
show that the radial (the strongest) component of the mag- 
netic field is very similar in both plots (see top left panel in 
Fig. 12 of D07). It shows a strong positive pole, sHghtly mis- 
placed relative to the rotational axis, then a negative octupo- 
lar ring above the equator, and a positive octupolar ring be- 
low the equator. In the D07 plot, the octupolar rings are not 
smooth in the azimuthal direction, which shows the presence 
of the non-axisymmetric modes. However, this difference is 
not very important for the present modeling, because most of 
the magnetic energy is in the axisymmetric multipoles. 

The surface distribution of the meridional component of 
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Model 


Ml 


Bi^ (G) 


Bo (G) 


po (g cm- 


-3) 


Mo (Moyr-i) 


Fo (erg cm ^) 


A^o (gcm^s-^) 


model 1 


1.5 


350 


5 


6.7 X 10- 


14 


4.7 X 10-9 


4.8 X 10^ 


2.8 X 10^6 


model 2 


1 


350 


7.5 


1.5 X 10- 


13 


1.0 X 10-s 


1.1 X 10^ 


6.2 X 10^^ 


model 3 


0.5 


350 


15 


6.0 X 10- 


13 


4.2 X 10-s 


4.4 X 10^ 


2.5 X 10^7 


model 4 


1.5 


700 


10 


2.7 X 10- 


13 


1.9 X 10-s 


1.9 X 10^ 


1.1 X 10^'^ 


model 5 


1.5 


1050 


15 


6.0 X 10- 


13 


4.2 X 10-s 


4.4 X 10^ 


2.5 X 10^^ 



Table 2. The reference values which depend on the parameters Bi^ and Jli. The dimensional values can be obtained by multiplying the dimen- 
sionless values obtained in different simulations by these reference values (see sec. 2 for details). In models 1-3, we varied the parameter Jli, and 
in models 4 and 5, the parameter Bi^. 




-1500 -750 750 1500 



Figure 2. The surface magnetic field in the dipole+octupole model 
of V2129 Oph = 1.5, 123 = 0.33, Bi = 30°, 82 = 20°, 03 

72' ) as seen from the equatorial plane (left panel), the north pole 
(middle panel) and the south pole (right panel) . The colors represent 
different polarities and strengths of the magnetic field; red/yellow 
means positive, blue/dark green means negative. 



/ii = 0.5Bi^Rl,we obtain the distance where the dipole and 
octupole components are equal: 

Here, w^e took into account that Bi^ = 0.35 kG and Bs^ — 1.2 
kG. That is, the octupolar field is expected to dominate at 
req < 1.6i?^, whereas the dipole field is expected to dominate 
at larger distances. 

Fig. [3] shows the distribution of the intrinsic magnetic 
field in our dipole+octupole model of V2129 Oph (at t = 0). 
The left panel shows that the field near the star has an oc- 
tupolar structure up to a radius of r ^ (1.5 — 2)i?^ which is 
in agreement with the Veq estimated above. The right panel 
shows that at larger distances the field has a more ordered, 
dipolar structure. 



the field. Be (Fig.[T] right panel) is quite similar to that in Fig. 
13 of D07 (left bottom panel). In both cases we see an inner 
ring of positive polarity (red), and an outer ring of negative 
polarity (blue) . The distribution of the azimuthal component, 
B(p, shown in Fig.[T] is somewhat similar to that shown in Fig. 
12 of D07 (left middle panel). However, our azimuthal field is 
weaker because we do not take into account the toroidal com- 
ponent of the field (see the right panel of Fig. 13 of D07 for 
the toroidal component) . We suggest that the toroidal and the 
higher-multipole poloidal components may influence the spot 
shapes near the star and can be included in future modeling. 

Fig. [2] shows the distribution of the modeled field on 
the surface of the star. The middle and right panels show 
strong northern (positive) and southern (negative) poles. The 
left panel shows that there are also northern (negative) and 
southern (positive) belts. One can see that the field is a typ- 
ical octupole field tilted by 20° (see also |Long7^omanova &, 
|Lamb|[2010a| ). The dipole component is present but is much 
smaller. 

An important radius is the radius where the dipole and 
octupole components of the field are approximately equal. Be- 
low, we derive an approximate value of such a radius, assum- 
ing that the dipole and octupole moments are aligned with 
the rotational axis. The magnetic fields of the dipole and oc- 
tupole vary with distance as (see L ong, Romanova & Lamb| 



2010a Eqn. 1): (i) in the equatorial plane: Bi 
B3 = l.S/is/r^, and (ii) in the polar direction: Bi 
B3 — A/iis/r^. Equating Bi and B3 in the equatorial plane 
and evaluating /ii and /is at the surface of the star (e.g.. 



MiA 

^ 2/ii/r 



4 MODELING OF ACCRETION ONTO V2129 OPH 

In simulating accretion onto CTTS V2129 Oph, we use its 
mass and radius derived from observations which are — 
I.35M0, — 2ARq (see D07 and references therein) and 
rotation period = 6.53 days (see Shevchenko & Herbst| 
1998 ) . These values imply that the corotation radius rcor — 
6.8R^, An estimate age of V2129 Oph is 2 x 10^ years (D07). 

A number of simulation runs were performed. To change 
the size of the magnetosphere (which is equivalent to the 
changing of the dimensional accretion rate at the fixed mag- 
netic field on the star), we varied the parameter Jli from 0.5 
to 5 and obtained a set of simulation runs where the disk 
stops at different distances from the star. We find from the 
simulations that at /Ii > 3 the inner disk radius is larger 
than the corotation radius and the disk moves away due to 
the "propeller" effect. At smaller parameters, 1.5 < /Ii < 3 
the magnetosphere is large, n > {6 - 8)R^, however, the ac- 
cretion rate is much smaller than the observed one. We chose 
some intermediate case with /Ii = 1.5 as the main one. In that 
case the disk stops far away from the star rt ~ 6i?^ (which 
is close to the corotation radius, Vcor ~ 6.8R^), and at the 
same time an accretion rate is not too low. We refer to this 
case as the main one. It corresponds to model 1 in Tab. [s] 
For this case we performed detailed analysis of different fea- 
tures of the flow, hot spots and the large-scale magnetic field 
evolution (see sections 4.1-4.3, 4.5, 4.6). In models 2 and 3 
(see Tab. [3| we investigate additional cases, where the disk 
comes closer to the star and accretion rate is higher and is 
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Figure 3. The distribution of the magnetic field in the dipole+octupole model of V2129 att = 0. The field is close to octupole near the star (left 
panel), but the dipole component dominates at r > 2R^ (right panel). The colors along the field lines represent different polarities and strengths 
of the field. 




Figure 4. 3D view of matter flow onto the star in the 
dipole+octupole model of V2129 Oph (model 1). The disc is shown 
by a constant density surface in green with p = 6 x 10~^^ gcm~^. 
The colors along the field lines represent different polarities and 
strengths of the field. The thick cyan, white and orange lines rep- 
resent the rotational axis and the dipole and octupole moments re- 
spectively. 



closer to the observed one. How^ever, the low-latitude octupo- 
lar belt spots become stronger than the polar spot (see sec. 
4.7). In models 4 and 5 w^e keep the same octupolar field, but 
increase the dipole magnetic field by a factor of 2 and 3 re- 
spectively. This helps to increase an accretion rate and at the 
same time octupolar belt spots are w^eak and the spot shape 
is closer to the observed one (see sec. 4.9). In model 6 we 
investigate accretion on to a star with a pure dipole field, and 
compare results with the main case (see sec. 4.4). Models 7-9 
are theoretical models (see sec. 4.8). 




p(10- g cm-0 



0,07 0.11 0.19 0.31 0.53 0.88 1.48 2.47 



Figure 5. Density distribution in different slices a.tt = 9 for model 1. 
Panels (a) and (b) show xz slices att = and t = 9. Panels (c) and 
(d) show yz and xy slices at t = 9. The red line shows the distance 
at which the matter stress equals the magnetic stress. 



4.1 Matter flow around the star 

Here we discuss results of simulations in our main case 
(model 1 in Tab. [3|. Initially, we place the inner edge of 
the disk at the radius r — 8R^. The small viscosity incor- 
porated into the disk leads to inward matter flow. First, the 
disk moves inwards over t ^ 4 rotations (Keplerian rotations 
at r = 1). Later, the funnel streams form and the matter ac- 
cretes onto the star quasi-stationary. We were able to reach 
time t = 8 — 10 in this type of simulations. Fig. [4| shows the 
matter flow and selected field lines at t = 9. One can see that 
the disk is stopped by the dipole component of the field, and 
matter flows towards the star from two directions, which cor- 
respond to the directions of the dipole component tilt (f^/xi 
plane) . 

Fig. [s] shows the density distribution in three slices, xz, 
yz and xy. One can see that the disk is disrupted by the 
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dipole component of the field far away from the star, at 
< rt < 7R^, and the funnel streams form at r 6R^. 
The red line shows the surface where the matter and magnetic 
field stresses balance each other. The dominant stress is given 
by the 0(/)-component of the stress tensor. The magnetospheric 
radius rm is then the radius where /3 = (p + pv^)/{B'^ /Stt) 



1. Inside this radius the magnetic stress dominates (|Ghosh &| 
Lamb|1978l ^l979 ; Long, Romanov a & Lambf 2010a; see also 
Bessolaz et al. 20 08 for different definitions of the inner disk 
radius). We obtain rm ~ (4 — 5)R^. 

Closer to the star, the octupole component of the field be- 
comes stronger, and each funnel stream splits into two parts, 
one of which flows towards the magnetic pole, and the other 
into an octupolar belt. This produces two types of accretion 
spots on the star. We discuss this in more detail in g4.3| 




p(10-^'gcm-^) 



0,67 0,84 1.04 1.30 1,63 2.03 2.54 3,17 

Figure 7. The calculated hot spots viewed from different directions 
Sitt = 9: from the equatorial plane (left-hand panel); from the north 
pole (middle panel); and from the south pole (right-hand panel). The 
color contours show the density distribution of the matter. Solid lines 
represent the magnetic moments of the dipole (/^i) and octupole 
(M3). 



4.2 Magnetic field structure and the validity of the 
potential approximation 

Next, we investigate the magnetic field evolution due to 
the disk-magnetosphere interaction. Fig. [6] (left two panels) 
shows the initial magnetic field distribution near the star (top 
panel) and in a larger region (bottom panel) . These plots cor- 
respond to the potential approximation which is often used 
for extrapolating the surface magnetic field to larger distances 
(e.g., D07). As the simulations proceed, the plasma and mag- 
netic field motions create currents and fields in the simulation 
region and the resulting magnetic field geometry strongly de- 
parts from the potential one. The dominant processes are con- 
nected with the differential rotation of the foot-points of the 
magnetospheric field lines which have different angular ve- 
locities at the star's surface, in the disk, and in the corona. 
Differential rotation leads to conversion of the rotational en- 
ergy into magnetic energy and to inflation of the field hnes 
(eXjAly 1980; Lovelace et al.|1995[ [Romanova et al.|1998| 
[Agapitou & PapaloizoU|2000| ) . The foot-points threading the 
star usually rotate faster than those threading the disc, and 
the magnetic field lines wrap around the rotational axis form- 
ing a magnetic tower (e.g. Lynden-Bell 1996) . Formation of 
such a tower has been observed in different numerical sim- 



ulations (e.g., Ustyugova et al.||2000 


Kato, Hayashi & Mat- 


sumoto|2004||Romanova et al.|2004b 


2009; Ustyugova et al. 


2006 ) and is a natural outcome of disk-magnetosphere inter- 



formation of such a tower. So, at large scales the magnetic 
field strongly departs from the potential one. 

The situation is different near the star. Fig. [6] (top panels) 
shows that the field structure at t = 9 is almost the same as at 
t = 0. This is the region r < (4 — 5)i?^, where the magneto- 
sphere is not disturbed by the disk, and magnetic stress dom- 
inates over the matter stress, and hence the magnetic fields 
produced by plasma flow are insignificant compared with the 
main magnetic field of the star produced by currents flowing 
inside the star. The potential approximation is valid in this 
region. 

We conclude that the potential approximation used in 
many works to recover the external magnetic field from the 
observed one (e.g. [Jardine et al.] [2002; D07) is sufficiently 
good in the vicinity of the star, where the magnetic stress 
dominates over the matter stress (inside the Alfven surface, 
r < rm), and where the magnetosphere is not disturbed by 



the accretion disc. However, the field strongly departs from 
the potential one at larger distances. 

It is often suggested that in young stars a strong stellar 
wind stretches the stellar field lines in the radial direction, 
like in the solar corona, and this determines the geometry 
of the external magnetic field (e.g. Safieryi998[ [Gregory et| 
|aL][2006[ [Jardin e, Greg ory & Donati[[2008| . We^should note 
that this solar-type scenario can dominate in the case of disk- 
less, weak-line T Tauri stars. However, if a large-scale mag- 
netic field truncates the accretion disk, then a different, mag- 
netic tower geometry dominates. Such a tower can acceler- 
ate a small amount of coronal matter up to high velocities 
and can carry some angular momentum from the star to the 
corona. However, both processes are significant only in cases 
of rapidly rotating stars (e.g., [Romanova et al.[[2005[ [Ustyu-[ 
[gova et a L] [2006j ), and are less important in slowly rotating 
stars, like V2129 Oph ( Romanova et al. 2009). 

G regory et al. ( 2008 ) noted that the dominance of the 
octupolar component near the star causes a smaller fraction 
of the stellar surface to be connected to open field lines com- 
pared with the case of a purely dipolar field. We observed that 
in V2129 Oph the amount of flux in open field lines depends 
on the position of the disk and on the strength of the dipole 
component which dominates at r > 2R^ and which deter- 
mines the amount of open flux. The field lines inside the inner 
disk are closed. All other field lines of the dipole component 
inflate and represent open magnetic flux. The octupolar com- 
ponent contributes towards open flux in polar region, but this 
flux is much smaller than the dipole open flux, as suggested 
by [Gregory et al.[ ( [2008] ). 



4.3 Accretion spots in V2129 Oph 

DO 7 analyzed the brightness enhancements observed in 
V2129 Oph in the chromospheric lines, like in Call IRT, and 
also in the Hel D3 line, which presumably forms in the shock 
wave near the star's surface. The right panel of Fig. 9 of D07 
shows the Call IRT excess emission relative to the basal chro- 
mospheric emission, which can be interpreted as the chro- 
mospheric "hot spot". Fig. 9 of D07 shows that the main 
part of the spot, the polar spot, is located at colatitudes 
0° < Oc < 50° with a center at Oc ^ 25°. The position of 
this spot almost coincides with the position of the magnetic 
pole (see Fig. 12 top left panel of D07). In addition, non- 
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Figure 6. Comparison of the initial (potential) magnetic field distribution att = (left panels) with the field distribution att = 9 (middle and 
right panels) . The top panels show the magnetic field in the vicinity of the star, while the bottom panels show the field distribution in the whole 
simulation region. The left and middle columns show the side view, while the right column shows the axial view of the field. 




Figure 8. The calculated hot spots in a polar projection down to 
120° from the north pole. Left panel: density distribution; right panel: 
energy flux distribution. The equator is shown as a solid line. The 
dashed lines represent the latitudes 30° and 60° respectively. 



axisymmetric spots of lower brightness are seen at higher co- 
latitudes up to Oc ^ 90°. 

In our numerical model, matter flows towards the star in 
funnel streams. We suggest that all matter which approaches 
the star freely moves inward through the star, and hence we 
neglect the complex physics of the stream-star interaction. In- 
stead, we calculate the matter and energy fluxes carried to- 
wards the star, and our accretion spot represents a slice in 
density/energy distribution across the funnel stream. 

In reality, the stream interacts with the dense regions of 
the chromosphere and photosphere, and forms a shock wave 
where gas is heated, and the energy is radiated mainly be- 
hind the shock front (e.g., "Lamzin"1995| |Calvet & Gullbring 
^98 ). The shock wave may oscillate due to cooling/heating 
instability behind the shock front ( [Koldoba et al.|20(58l ), and 
the stream-star interaction may lead to formation of waves on 
the star's surface and possible outflows ( [Cranmer |2008| [2009{ 
see also |Brickhouse et al.|2010l ). We suggest that the heating 



Y 




Figure 9. The color background shows the distribution of the total 
magnetic field on the surface of the star. The black contours show 
the energy distribution in the accretion flow at the star. The black 
solid lines show the direction of the dipole and octupole moments. 
The dashed line shows the direction towards the brightest part of 
the spot. The vertical dashed line shows the origin line for the phase 
count in D07 paper (clockwise). 



processes in the dense layers of star's atmosphere may excite 
such lines as Call IRT, Hel D3 and others, and hence we ex- 
pect a correlation between the position of the chromospheric 
spots observed by D07, and accretion spots observed in our 
simulations. 

Fig.[7|shows different views of the density spots. One can 
see that there are two main spots, the polar ones, which al- 
most coincide with the octupolar magnetic poles. There are 
also two elongated spots at much lower latitudes where mat- 
ter flows towards the magnetic octupolar belts. 
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Fig. [s] shows a polar view of the spots. The left panel 
shows the density distribution. The right panel shows the 
distribution of the energy flux which is calculated as: F = 
pv-r[{v'^ —vl)/2-\-^p/{'j — l)p] (here^;^ is the surface velocity 
of the star) and shows the location of the main energy depo- 
sition in the funnel stream. One can see that the polar density 
spot is located at colatitudes of 5° < < 45° and is centered 
at Oc ~ 30°. The octupolar belt spots are azimuthally elon- 
gated and are located at 60° < < 80°. The energy spots 
have similar shapes and locations. We also calculated the dis- 
tribution of the matter flux pVp and temperature T (calculated 
in the suggestion that all kinetic energy is converted into black 
body radiation) . We observed, that their distribution is similar 
to that of the density and energy fluxes shown in Fig.js] Any of 
these values can be used for analysis of accretion spots on the 
star. Note that the size of spots in our simulations depends on 
the chosen density/energy level, so that the densest/brightest 
parts of spots have a smaller size. 

Comparison of the accretion spots obtained in our sim- 
ulations with the Call IRT chromospheric "spots" observed 
in V2129 Oph (Fig. 9 of D07) shows that the polar spot is 
present in both cases and has a similar location and shape. 
The simulated spot is somewhat more compact than the ob- 
served one (see however §3.6 where we discuss the depen- 
dence of the spot size on the density/energy levels) . Also, the 
center of the observed spot is located at a slightly higher lat- 
itude than the simulated spot. In both cases, the polar spots 
are located near the magnetic poles seen in Fig. 1 of our pa- 
per (and Fig. 12 of D07), though the position of the simulated 
density spot is somewhat different in phase compared with 
the magnetic pole position. 

Fig. [9] shows that in our simulations an accretion spot 
(shown in black contour lines) does not coincide with the 
maximum of the magnetic field (red region on the star). A 
spot is leading (about 0.15 in phase). Here we note that the 
spot is located between the dipole and octupolar poles. We 
suggest that the main accretion stream moves toward the 
dipole magnetic pole and hence the accretion spot is located 
close to this pole. However, an octupole field also influences 
the accretion stream, and as a result a spot is located between 
the dipole and octupolar magnetic poles. This may explain 
why we observe the phase difference of about 0.15 between 
the main (octupolar) magnetic pole and accretion spot. The 
uncertainty in the tilt and phase of the dipole field derived 
from observations is about 10%. If we assume that the phase 
difference between the dipole and octupolar poles is smaller, 
then we expect that the accretion spot will be closer to the oc- 
tupolar magnetic pole (as observed by D07). There is also an- 
other possibility, which is that the funnel stream is dragged by 
the rapidly rotating disk and moves slightly ahead of the more 
slowly rotating magnetosphere, and this leads to the phase 
difference between the spot position and the magnetic pole 
(as shown in simulations by [Romanova et al.|2003|[2b04a| ). 

4.4 Comparison with a pure dipole case (model 6) 

We performed a special set of runs where the magnetic field 
of V2129 Oph is approximated with the dipole component 
only (zero octupolar field). The strength of the dipole field, 
0.35 kC, and its tilt, 6 = 30°, are the same as in the main 
dipole + octupole case. Figure [To] (top right) shows the matter 
flow and spots in the x2;-plane in the case of a pure dipole 



field. One can see that matter flows in two narrow streams 
and the hot spots have the azimuthally elongated, crescent 
shapes, typical for stars with the dipole field ( ||Romanova et al^ 
2004a; Kulkarni & Romanova 2005) . The spots are centered 
at a much lower latitude, Os ~ 45°, compared to spots in the 
dipole octupole case. 

The left panels show the xz cross-section and spots in 
our main dipole + octupole case for comparison. One can see 
that initially matter flows along the dipole field lines. How- 
ever, closer to the star, the octupole component determines 
the flow pattern: the funnel stream becomes wider and splits 
into a polar stream and an equatorial belt stream. Note that 
the polar stream is redirected by the octupolar component to 
much higher latitudes compared with the stream in a pure 
dipole case. 

We conclude that the high-latitude round accretion spots 
obtained in simulations of the dipole + octupole model are 
in much better qualitative agreement with the observed ac- 
cretion spots on V2129 Oph (D07) than the lower-latitude 
crescent-shaped accretion spots obtained in the pure-dipole 
model. 

In this work we have neglected the higher-order multi- 
poles which can dominate near the stellar surface. We should 
note that they cannot change the position of the main accre- 
tion spot, but instead can only redistribute the density inside 
the main spot. 



4.5 Area covered by hot spots 

Here we calculate the area covered by the accretion spots. We 
choose a density p, and calculate the area A{p) covered by the 
part of the spots with density higher than p, and the fraction 
of the star covered by these spots: /(/>) = A{p)/4:7iRl. We 
calculate this fraction for different values of p. To separate 
the contribution from the polar spots, we calculate separately 
the fraction of the star covered by all spots (ftot), and by po- 
lar spots only, fpoiar- We also calculate the fraction of the star 
covered by spots in the pure dipole case, fdipoie- Fig. 11 shows 



how the spot coverage varies with p. One can see that the area 
covered by spots in the dipole case is larger than that in the 
dipole + octupole case if p > 1.7 x 10~^^g cm~^ (p > 0.26). 
The situation is reversed if the density level is smaller. Note 
that in earlier comparisons of accretion to stars with dipole 
and quadrupole fields, we observed a similar trend: at higher 
density levels, the spots in the pure dipole case are larger than 
those in the dipole plus quadrupole cases, while at lower den- 
sity levels, the situation is reversed ( [Long, Romanova & Lamb| 
|2010a| ). 

In another figure. Fig. [12] we choose several density lev- 
els and show the size of spots corresponding to these den- 
sity levels. One can see that at the lowest density level, 
p = 7 X 10~^^g cm~^ (p = 0.1), the spots occupy ftot ^ 20%, 
while at the highest density level, p — 2.1 x 10~^^g cm~^ 
(p = 0.3) , they occupy ftot ~ 3%. If only the polar spot is 
taken, then the coverage for the same density levels is smaller 
— fpoiar ~ 8.0% and fpoiar ~ 2.3% respectively Fig. [TT 
shows that at even higher density levels, the spot coverage is 
very small. For example, at p = 3 x 10~^^g cm~^, ftot ~ 1%. 

Thus, the spot size strongly depends on the density level. 
This density distribution shows the structure of the funnel 
streams: the density is higher in the middle of the stream, 
and decreases toward the periphery. The matter flux, kinetic 
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Figure 10. Comparison of matter flow and hot spots in the dipole+octupole model of V2129 Oph (left panels) with a pure dipole case (right 
panels) at t = 9. The top panels show the density distribution in the xz— plane (color background) and sample field lines. The bottom panels 
show the density distribution in the hot spots as seen from along the y— axis (panels a,c) and along the rotational axis (panels b,d). 



energy flux and temperature have a similar spot-centered dis- 
tribution and hence the spot coverage is different at different 
energy/temperature levels with the hottest spots occupying a 
smaller fraction of the star's surface (see also |Romanova et al] 
|2004a| ). 

The size of the accretion spots reconstructed from obser- 
vations, e.g., of the Call line, also depends on brightness as 
seen in Figs. 9 and 11 of D07 and Fig. 6 of DIO. The fraction 
of the star covered with each of visible spots is approximately 
5% which is about 10% for both spots (observations take into 
account only the spots in the visible hemisphere, and do not 
take into account a probable antipodal spots; to compare sim- 
ulations with observations, we multiply the observed fraction 
by a factor of two). One can see that the brightest parts of 
the observed spots cover a much smaller area of the star. The 
position and shapes of the observed and simulated spots are 
similar: both are round and are located at the co-latitude of 
Oc ^ 30°. 




p(10-''g cm-') 



Figure 11. Fraction of the star's surface covered by spots of density 
p and higher att = 9. The bold line shows the spot coverage in the 
dipole -I- octupole model of V2129 Oph where both polar and octupo- 
lar belt spots are present. The thin solid line shows same, but for the 
polar spot only. The dashed line shows same but for a pure dipole 
model ofV2129 Oph. 



Fig. [12] (middle panel) shows that at this spot coverage 
(/^o^ = 10%) the polar part of the spot strongly dominates, 
while the octupolar belt spots are weaker. We conclude that 
there is a good correspondence between the sizes and lon- 
gitudinal positions of the main, brightest parts of the spots 
obtained in observations and simulations. However, there is a 
difference in shapes and phases of the lower-density/energy 
parts of the spots (see also §4.4 and §4.6). 

It was predicted that the area covered by spots in the 
case of complex magnetic fields should be smaller compared 
with the spot area in the pure dipole case (e.g., Mohanty & 
Shu 2008; Grego ry et al.||2008| see also Calvet & Gullbring 
1998) . In addition, [Calvet & Gullbring, ( ,199 8) derived from 
their model and observations of a number of CTTSs that their 
spot coverage is expected to be very small, less than a per- 
cent in most stars. However, our simulations show that the 
spot coverage is higher than that predicted by [Calvet & Gull-| 
[bringl ( [1998| ) . It is possible that these authors analyzed only 
the central, brightest parts of the spots. This problem requires 
further investigation. 




0.6 OA 
Rotational cycle 

Figure 13. The rotationally modulated longitudinal magnetic field 
derived by DO 7 from observations of the emission line Call IRT (blue 
line) compared with the rotationally modulated radiation from spots 
obtained in our simulations (red line) and from simulations of the 
pure dipole field (dashed black line) for inclination angle i = 45°. 
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Figure 12. Example of spots with a cut at the different density levels p and corresponding spot coverage for the full spot ftot- 
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Table 3. The Table combines different observational properties of V2129 Oph (top three rows), and shows different models used in simulations 
and their results (see text). The dimensionless parameters Jli and JI3 determine the dimensionless properties of the magnetosphere; Bi^(G) 

and Bs^ (G) are the polar dipole and octupole fields given by observations; M is the dimensionless parameter of the accretion rate derived from 
simulations; M(M0yr-i) is the dimensional accretion rate; rt/Ri, is the truncation radius of the disk given in units of stellar radius. The last 
column describes the properties of spots at such a density/energy level when they cover about / = 10% of the stellar surface. The top three 
rows show the observed accretion rates including its maximum, minimum and middle values. Model 1 is the "main" (reference) model and it is 
described in detail in different subsections of sec. 4. It is also used for comparisons with other models. Models 2 and 3 show results of simulations 
in cases where accretion rate is higher, and the disk comes closer to the star (see sec. |4.7| > . Models 4 and 5 show results of simulations in cases, 
where the dipole component is twice and thrice as large compared with observations (see sec. |4.8| >. Model 6 shows result of simulations in a pure 
dipole case (see sec. |4.4| >. Models 7-9 describe results of the theoretical model (see sec. |4.7| >. 



4.6 Rotational modulation 

DO 7 observed that the fluxes and amplitudes of the Zeeman 
signatures in the Call IRT and Hel D3 emission lines exhibited 
strong rotational modulation over the rotational cycle of the 
star w^hich is close to the 6.53 d period in V2129 Oph. The 
time-dependent longitudinal field derived from these emis- 
sion lines also show^s the rotational modulation (see Fig. 5 of 
D07). The emission in Call IRT traces the regions of higher 
excitation of this line in the chromosphere, while the Hel D3 
line is thought to form in the shock w^ave ( [Beristain, Edwards] 
|& Kwan|1998| ) . These features can be compared with the ro- 
tationally modulated accretion spots obtained in our simula- 
tions. 

We calculate the rotationally modulated radiation from 
accretion spots towards the observer located at an angle i rel- 
ative to the rotational axis of the star. To calculate the radia- 
tion, we suggest that all the energy flux of the funnel stream 
is converted into isotropic black-body radiation in the shock 
wave (see details in |Romanova et al.|2004a| ) . We choose some 
late moment of time (t = 9) when the position of the spots 



does not vary much and obtain the rotationally modulated 
radiation from such spots. 

Fig. [13] compares the rotationally modulated longitudi- 
nal magnetic field derived from Call IRT observations by D07 
(blue line) with our rotationally modulated spot curves (red 
hnes) which we calculated at the inclination angle i = 45° 
(D07). One can see the gradual rise and sharper drop in both 
simulated and observed curves which may be a sign of the 
spots asymmetry. We calculated a similar rotationally modu- 
lated curve for V2129 Oph with the magnetic field approxi- 
mated by a pure dipole field and found that the curve is of 
a sinusoidal shape (see the dashed black line in Fig. [Ts] ) . We 
suggest that the non-sinusoidal shape may be connected with 
the presence of the low-latitude (octupolar belt) part of the 
spot in our simulations, and with the low-latitude components 
of the spots in the observations. 

The observed and simulated curves have a phase differ- 
ence of about 0.1-0.2 which probably reflects the phase dif- 
ference between the simulated and observed spots. 
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Figure 14. Result of simulations at different magnetospheric sizes 
(parameters Jli) for models 1, 2, 3 (see Tab. |3}. Top panels show 
density distribution (background), selected field lines (yellow lines), 
and /3 = 1 line (red line) in the xz-slices. Bottom panels show the 
energy flux distribution on the star in the polar projection. 



4.7 Accretion rate onto V2129 Oph 

The accretion rate in V2129 Oph has been estimated by D07 
using the Call IRT l ine to be M 4 x 1Q~^ Mq yr~^ (us - 
ing the approach of |Mohanty, Jayawardhana & Basri||2005| ). 
Using fluxes from several emission lines and more recent em- 
pirical relations between line emission fluxes and accretion 
luminosities ( [Fang et al.|2009| ), a revised mass accretion rate 
was derived from the same spectra: \og{M / {Mq yr~^)) = 
-9.2 ± 0.3, or: M = 6.3 x 10~^^ Mq yr"^ with the maxi- 
mum and minimum values of Mmax — 1.2 x lO^^Moyr"^ 
and Mmin = 3.2 x lO~^°M0yr"\ 

Below we compare the observed accretion rate with the 
rates obtained from simulations and derived from theoretical 
estimations. The dimensional accretion rate is 

AW., . 5Ai,.l.lxI0-M^(|)(^)' 

( ^* VV^V (5) 

In the case of /Ii = 1.5, we take M ^ 0.015 from the simula- 
tions, and Mo from Tab.0and obtain a dimensional accretion 
rate Msimui = 7.2 x lO^^^l^oyr"^ . This accretion rate is about 
9 times smaller than the observed accretion rate (see Tab.js]). 

In an attempt to resolve this issue, we performed simula- 
tion runs at smaller values of the magnetospheric parameter, 
Jii — I and 0.5 (see models 2 and 3 in Tab.[3|. In each case, 
we adjusted the dimensionless octupolar moment using the 
ratio: Ji^/Jli — 0.22, thus fixing the ratio between the dipole 
and octupole moments. We observed from the simulations 
that for /IT = 1, the accretion rate is higher and it matches the 

minimum value derived from observations, Mmin (see model 
2 in the Table). For Jii — 0.5, the accretion rate is sufficiently 
high to match the observed accretion rate (see model 5) . 
Fig. [l4| (top panels) shows matter flow at different Jn . 



One can see that the size of closed magnetosphere systemat- 
ically decreases with /II. We determine (very approximately) 
the truncation radius as the radius at which the last closed 
(or, strongly disturbed) field line crosses the equatorial plane 
(shown as a dashed line on the plot). Then for /Ii = 1.5, 1, 0.5, 
we obtain a set of radii, rt ~ 6.2, 4.3, 3.4. Each of these radii 
has an error of about dz(0.1 — 0.3), because there are usually 
a few field lines which can be considered as slightly disturbed 
lines. Nevertheless, we see a clear trend of rt decreasing with 
Ml. 

Fig. [14] (bottom panels) shows that an area covered with 
the octupolar belt spots increases when Jn decreases, because 
the disk is truncated closer to the star, where the octupolar 
component dominates. Note that in the case of /Ii = 1, the 
polar spot dominates over octupolar ring spots if we take the 
brightest/densest 10% of the stellar surface. However, in the 
case of /ii = 0.5, octupolar ring spots give the brightest part at 
any density/energy level. Such octupolar spots have not been 
observed by D07. So, the intermediate case (/Ii = 1, model 
2) can possibly explain the observations, though the trunca- 
tion radius rt 4.3 is somewhat smaller than the corotation 
radius, Rcor = 6.8i?^. 

To further investigate this issue, we estimated the accre- 
tion rate using a theoretical approach. Prior simulations done 
in our group show that the formula for the Alfven radius ( [El-| 
Isner &Lamb|1977| ): 

Tt = A;(G'M.)"'/'M"'/Vt^'. (6) 

is approximately applicable for the truncation of the accretion 
disks, lik 0. 5 ([Long, Romanov a & Love lace|2005 | see also 
[Bessolaz et aL||2008| for different approaches for finding the 
inner disk radius) . We can rewrite this formula in the conve- 
nient form: 

V1.35M0; \2ARq) yr 

Tab. [3] (models 7-9) shows theoretically derived values of M 
for the truncation radii n obtained in models 1-3. One can see 
that these accretion rates approximately correspond to those 
obtained in the 3D simulations. 

Hence, both simulations and theory point to the fact that 
a star with an equatorial dipole magnetic field of 175 G can 
truncate the disk at the corotation radius only if the accretion 
rate is very small, smaller than that found from observations. 
Therefore, the disk should be truncated at smaller distances 
from the star to explain the accretion rate. 

4.8 Test simulation runs at larger dipole field 
component (models 4 and 5) 

To further resolve this issue, we suggest that the dipole com- 
ponent is larger than that reconstructed by D07 from obser- 
vations. In fact, the July 2009 observations revealed a dipole 
component about three times stronger (DIG). So, our test 
simulations may be relevant to some other moments of time, 
when the dipole field is stronger. 

We performed two test simulation runs (models 4 and 
5) where the octupolar component is fixed at B^^ — 1.2 
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Figure 15. Test simulations of cases where the ratio of the dipole 
to octupole component is twice and thrice as large compared with 
the main case (from left to right: models 1, 4 and 5; see also Tab. 
|3}. Top panels show density distribution (background), selected field 
hues (yellow lines), and ^ = 1 line (red line) in the xz-slices. Bottom 
panels show the energy flux distribution on the star in the plane pro- 
jection. Note that in each case we chose such energy level, that the 
spots coverage is ftot ~ 10% and hence it is close to one observed 
in V2129 Oph. One can see that the polar spots dominate in both 
models 4 and 5. 



kG, v\^hile the dipole component is tw^ice and thrice as large: 

Bi^ = 700 G (/Is/mi = 0.165) and Bi^ = 1050 G (/I3/M1 = 
0.11). 

Fig. [15] shows the results of these runs versus the main 
case. One can see that in all cases the disk stops at approx- 
imately same distance from the star (because Jii = 1.5 in 
all cases) . However, the dimensional accretion rate is propor- 
tional to Bf^ and hence it is expected to have about 4 and 
9 times higher accretion rate compared with the main case. 

Exact value also depends on the dimensionless parameter M 
obtained from simulations. Tab. [s] (row 12) shows that the 
accretion rate in these two cases is sufficiently high. One can 
see that even the dipole field of Bi^ — 700G is sufficient in 
disrupting the disk at far distances from the star, and at the 
same time to give the observed accretion rate. Another ad- 
vantage of simulations with a higher dipole component is that 
the main polar spot strongly dominates which gives a better 
match with observation of accretion spots in V2129 Oph. 

4.9 Calculation of the torque and spin evolution 

We also calculate the torque on the star from the disk- 
magnetosphere interaction. Matter coming from the disk in 
funnel streams can move faster or slower than the star and 
hence it creates either a positive or negative azimuthal field 
on the star, and correspondingly a positive (spin-up) or neg- 
ative (spin-down) magnetic torque, Nf. In addition, matter 
brings some positive angular momentum which spins the star 
up, Nm- We calculated both the magnetic, Nf, and matter, 
Nm, torques, and found that the torque associated with the 
magnetic field is about 100 times larger than that associated 
with matter. This is typical in accretion through the funnel 



streams (see also [Romanova et aLl|2002| |2003| ). We calcu- 
lated the torque for models 1-6 and the results are presented 
in Table [4] We observed from the simulations that the torque 
is small and it wanders around zero taking both positive and 
negative values in models 1, 5 and 6 (see Tab.|4]). In models 
2-4, the torque is slightly positive and the^star spins up. We 
calculated the dimensional torque Nf — NfNo using corre- 
sponding reference values A^o from Tab. [2] 

Now, we can estimate the time-scale of spin evolution. 
The angular velocity of V2129 Oph is (at its period of P = 
6.53 days) Q = 27t/P ^ l.lx 10~^s~^, its angular momentum 
is J = kM^RlQ = 8.5 X lO^^k gcm^s, where k < 1, Taking 
A: = 0.4 and different values of Nf for different models from 
Tab. [4] we obtain the characteristic time-scale as r = J/Nf 
(see Tab. [4]). One can see that the time-scale in all models 
is longer than the age of V2129 Oph which is estimated to 
be 2 X 10^ years (see D07). Hence, the torque observed in 
simulations is not sufficient to spin down the star to its present 
slow rotation. 

If the star has lost it spin in the past due to disk- 
magnetosphere interaction, then we should suggest that in 
the past the accretion rate was higher (to ensure rapid spin- 
down), and the star's magnetic field was stronger (to ensure 
disk truncation at a large distance from the star). In fact it 
is expected that in the past, V2129 Oph has been fully con- 
vective and its magnetic field has been amplified due to the 
dynamo mechanism. The "propeller" mechanism can be re- 
sponsible for the efficient spinning down at this stage (e.g., 
[Romanova et al.|[2005t ?). The stellar winds can also carry 
angular momentum out at different stages of star's evolution 
(e.g., [Matt & Pudritz|200"5l [20081 [tonmerf 2008'/2^^ 

On the other hand, the July 2009 observations by DIO 
did show that the magnetic field of the star is stronger than 
in 2005 observations, and maybe the field of V2129 Oph is 
higher on an average than in the June 2005 observations. 
Note, however, that the torque is still somewhat small, even 
when the dipole field is three times stronger (see Tab. [4| for 
model 5). 

Note that if a torque at the present epoch is small, and 
a star is not in the rotational equilibrium, then the trunca- 
tion radius can be at any place (not at the corotation radius) . 
In such a case one of the main condition which should be 
satisfied in comparisons of models with observations is the 
similarity of the spot shapes and positions. 



5 CONCLUSIONS AND DISCUSSIONS 

We performed global 3D simulations of disk accretion onto a 
star with parameters close those of T Tauri star V2129 Oph. 
The large-scale magnetic field of V2129 Oph is approximated 
by slightly tilted dipole and octupole magnetic moments with 
polar fields of 0.35 kG and 1.2 kG respectively. Test simu- 
lations for the case of a pure dipole field, smaller magneto- 
spheres and a larger dipolar field component were performed 
for comparison. Below, we summarize our findings and com- 
pare our results with the D07 observations/model: 

1. The field lines connecting the disk and the star inflate and 
form magnetic towers above and below the disk. The poten- 
tial (vacuum) approximation is still valid inside the Alfven 
(magnetospheric) radius, rm < (4 — 5)i?^, where the mag- 
netic stress dominates over the matter stress. At larger dis- 
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Model ofV2129 Oph 


Ml 


M3 




Nf (g cm^s-^) 


r (yr) 




model 1. dip+oct 3D 


1.5 


0.33 


±0.006 


±1.7 X 10^4 


6.5 X 10^ 


rot. equilibrium 


model 2. dip+oct 3D 


1.0 


0.22 


+ 0.01 


±6.2 X 1034 


1.7 X 10^ 


spin-up 


model 3. dip+oct 3D 


0.5 


0.11 


+ 0.015 


±3.7 X 1035 


2.9 X 10'^ 


spin-up 


model 4. dip+oct 3D 


1.5 


0.165 


+ 0.01 


±1.1 X lO^s 


9.7 X 10'^ 


spin-up 


model 5. dip+oct 3D 


1.5 


0.11 


±0.02 


±5.0 X 10^^ 


2.1 X 10'^ 


rot. equilibrium 


model 6. dip 3D 


1.5 


0.00 


±0.003 


±8.3 X 1033 


1.3 X 10^ 


rot. equilibrium 



Table 4. The Table shows the dimensionless torque Nf, dimensional torque Nf = NoNf, and the time-scale of spinning-up/down r for models 



tances, the magnetic field distribution strongly departs from 
the potential one (see Fig.|6|. 

2. Simulations show that the disk is disrupted by the dipole 
component of the field and matter flows towards the star in 
two funnel streams. The flow in the stream is redirected by 
the octupole component with part of matter flowing towards 
the high-latitude pole, and another part into the low-latitude 
octupolar belt. The polar spots are located at high latitudes 
even if the disk is truncated close to the star (say, at r ^ 

0.5i?cor). 

3. In the dipole case the spots are latitudinally-elongated, they 
have a typical crescent shapes and are located at lower co- 
latitudes Oc ~ 45° (see also R04). The polar chromospheric 
spot observed in V2129 Oph (DO 7) is much closer in shape 
and position to the round, high-latitude polar spot obtained 
in the dipole + octupole model than the crescent-shaped spot 
obtained in the pure dipole model. 

4. If the disk is disrupted at rt ~ 6.2R^, like in our main 
case (model 1), then the accretion rate obtained from the 
simulations is about an order of magnitude smaller than 
that obtained from observations. If the disk is disrupted at 
rt ~ 3AR^ (model 3), then the accretion rate has a good 
match, but the low-latitude octupolar belt spots strongly dom- 
inate, which is in contrast with observations. However, if the 
disk is disrupted at a somewhat larger distance from the star, 
rt ~ A.SR^ (model 2), then the brightest part of the spot 
is located in the polar region, while the octupolar ring spots 
are much weaker. In this case the accretion rate is smaller 
than the observed accretion rate, but it is within the error bar 
given by DIO (see Tab. [3|. Hence, model 2 can possibly ex- 
plain observations of V2129 Oph. In that case, future, more 
sensitive observations should reveal the arcs of the octupolar 
belt spots. 

5. Both simulations and theoretical estimations point to the 
fact that a weak dipole field of Bi^ = 350G cannot disrupt 
the disk at large distances from the star, unless the accretion 
rate is very low In a test simulation with a dipole field twice 
as strong (Bi^ = 700 G, model 4), the accretion rate matches 
the observed accretion rate and the polar spots dominate. This 
model explains observations better than the other models. 

6. We observed that the main (polar) accretion spot is shifted 
relative to the main (octupolar) magnetic pole (at the phase 
of 0.15) towards the dipole magnetic pole (see fig. |9|. Such 
a shift may be connected with the fact that matter flows to- 
ward the dipole magnetic pole and is only partially redirected 
towards the octupolar pole. The dipole and octupole poles dif- 
fer in phases, and this may lead to the observed phase differ- 



ence. On the other hand, the funnel stream can be shifted for- 
ward/backward relative to the stellar surface and the phase 
of the accretion spot may differ from that of the magnetic pole 
( [Romanova et aL|2003]|2004a| ). 

7. The torque on the star is small and corresponds to a 
small accretion rate. The time-scale of spin evolution (r ^ 
2 X 10^ — 10^ years) is longer than the estimated age of V2129 
Oph of 2 X 10^ years. Hence, disk-magnetosphere interaction 
cannot be responsible for the star's spin evolution at present 
epoch. Hence, we conclude that the star probably lost most 
of its angular momentum during the early stages of its evolu- 
tion, for example, at the stage when a star has been fully con- 
vective and had a stronger, dynamo-generated magnetic field. 
The "propeller" mechanism could be responsible for the rapid 
spin-down ( [Romanova et aL|2005t|Ust3aigova et al.|2006| ). In 
addition, at all stages of evolution angular momentum flows 
into stellar winds (e.g. [Matt & Pudritz|2005l|2008| ). 

Among the different simulation runs presented here we 
prefer those runs where the disk is disrupted at a large dis- 
tance from the star, within (0.7 — 0.9)i?cor, in which the star is 
expected to be in the rotational equilibrium state (zero torque 
on average) (e.g., ,Long, Romanova & Lovelace]|2005 j . Su ch 
a state is probable in slowly rotating CTTSs ( Konigl|1991[ ?; 
'Romanova et aL|200 2) . However, the low values of torque ob- 
tained in our simulations raise the question of whether V2 129 
Oph is in the rotational equilibrium state, and hence the po- 
sition of the inner disk is less important for the analysis. 

A new set of observations of V2129 Oph done in July 
2009 have shown that the octupole and dipole components 
are about 1.5 and 3 times larger compared with the June 
2005 observations discussed in this paper (DIO). If the field 
of V2129 Oph varies due to the dynamo mechanism, then (at 
a steady accretion rate) the disk will move closer to the star 
during periods when the field is weaker, and further away 
when it is stronger. Then, model 2 may explain the June 2005 
observations, and models 4 or 5 (with dipole moments twice 
and thrice as large) the July 2009 observations. If the dipole 
component is large on an average, say, about (0.7-l)kG, then 
the accretion is sufficiently high to explain observations. 

[Jardine, Gregory & Donati| ( [2008| ) investigated the prop- 
erties of the funnel streams in V2129 Oph for a fixed magnetic 
field geometry. They chose a much higher accretion rate in 
their model M ^ lO^^Mo/yr, which is a typical accretion rate 
for mildly accreting CTTSs (see also D07, [Eisner et al.|2005) . 
The density they obtained in the funnel streams was about a 
hundred times higher compared with our model. At this ac- 
cretion rate, however, the disk will be truncated very close to 
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the star unless the dipole magnetic field is even higher than 
discussed above. 

There is an interesting possibility that most of the disk 
matter may accrete directly through the equatorial plane due 
to magnetic interchange instabilities ( Romanova, Kulkarni &| 
|Lovelace| [20081 [Kulkarni & Romanova] |2008| |200 9L while 
a smaller amount of matter climbs the magnetosphere and 
forms accretion spots. This model can possibly explain the 
difference between the low accretion rates obtained by DIO 
from spectral lines, and the much higher accretion rates de- 
rived from the disk measurements (e.g., [Eisner et al.|2005j ). 

We performed similar 3D MHD simulations of accretion 
onto another T Tauri star, BP Tau ( Long et al.|2010b| ) where 
the dipole and octupole components of 1.2 kG and 1.6 kG 
strongly dominate over other multipoles ( [Donati et al.|2008| ). 
In that case the dipole component is much stronger than in 
V2129 Oph, and there is a better match between the sim- 
ulated and observed accretion rates. Both V2129 Oph and 
BP Tau are examples of stars with complex fields where the 
dipole component determines the disk-magnetosphere inter- 
action, while octupolar component shapes the spots. 

In some T Tauri stars the magnetic field is more complex 
than inV21 29 Oph, like in CTTSs CVCha and CRCha ( |Hussain[ 
jet al.|2009l ). Modeling of such stars requires incorporation of 
even higher order multipoles, which can be done in future 
research. 
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